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ABSTRACT 


The  objective  of  this  research  is  to  model  and  study  the  effects  of  transbasin 
internal  waves  on  low-frequency  signal  transmission  through  the  South  China 
Sea  (SCS)  basin.  Specifically,  the  fluctuations  in  the  multipath  arrival  structure  of 
a  400-Hz  acoustic  pulse  transmitted  through  a  distance  of  167-km  in  the  SCS 
basin  in  the  presence  of  an  internal  ocean  soliton  was  modeled  and  examined. 
The  modeling  entailed  the  integration  of  a  raytracing  program  with  an  eigenray 
search  and  arrival-structure  calculation  program,  and  the  use  of  measured 
bathymetry  and  inferred  bottom-loss  characteristics  from  previous  research.  A 
range-dependent  perturbation  was  added  to  a  range-independent  background 
sound-speed  profile  to  model  the  varying  sound-speed  field  as  the  nonlinear 
ocean  internal  soliton  propagates  along  the  transmission  path.  All  cases  studied, 
each  simulating  a  soliton  at  a  different  location,  had  six  distinct  acoustic  arrivals 
that  suffered  from  large-amplitude  fluctuations  (~  10  dB).  The  factors  that  affect 
the  amplitude  of  the  arrivals  are  changes:  in  the  number  of  bottom  interaction,  in 
ray  tube  spreading,  phase  interference  and  in  the  number  of  eigenrays  making 
up  an  arrival.  The  results  also  show  that  the  closer  the  soliton  to  the  receiver,  the 
less  impact  the  soliton  has  on  the  arrival  structure. 
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I.  INTRODUCTION 


A.  NAVAL  IMPORTANCE 

The  current  Defense  Strategic  Guidance  states  while  the  U.S.  military  will 
continue  to  contribute  to  security  globally,  we  will  of  necessity  rebalance  toward 
the  Asia-Pacific  region.  In  order  to  do  this,  the  Navy  must  be  able  to  understand 
the  operational  environment  within  the  Asia-Pacific  region.  The  South  China  Sea 
(SCS)  is  an  area  of  significant  strategic  and  tactical  importance  to  the  Navy.  A 
large  portion  of  the  world’s  commerce  travels  on  the  shipping  lanes  located  in  the 
SCS.  An  accurate  model  of  the  phenomenology  of  sound  propagation  through 
the  SCS  needs  to  be  developed  to  help  maintain  tactical  superiority.  This  model 
will  provide  valuable  information  for  operational  planning  and  execution  in  the 
South  China  Sea. 

B.  SOUTH  CHINA  SEA  BACKGROUND 

The  SCS  is  the  westernmost  marginal  sea  of  the  Pacific  Ocean.  The  SCS 
has  borders  of  Straits  of  Malacca  to  southwest,  the  Taiwan  Strait  to  the 
northeast,  and  the  Luzon  Strait  to  the  east.  The  SCS  stretches  from  Taiwan  to 
Singapore.  The  bottom  consists  of  a  deep  basin  with  depths  greater  than  4,000 
m  and  two  extended  continental  shelves,  one  to  the  north  and  one  to  the  south. 

The  SCS  has  environmental  variability  on  numerous  scales  ranging  from 
super-tidal  to  inter-annual.  These  include  intrusion  of  the  Kuroshio  Current, 
variable  wind  forcing,  and  large  amplitude  internal  waves.  The  seasonal 
variability  is  governed  by  the  monsoons.  The  variable  wind  forcing  is  due  to  the 
tropical  monsoon  climate.  The  Kuroshio  Current  is  the  western  boundary  current 
of  the  northwestern  Pacific  Ocean.  It  has  been  observed  to  enter  via  the  Luzon 
Strait  and  then  have  widely  variable  behavior.  Possible  scenarios  include  simply 
exiting,  bifurcation,  or  forming  a  warm  mesoscale  eddy  that  can  propagate 
across  the  SCS  (Ramp  et  al.  2004).  The  large  amplitude  internal  waves  are 
generated  by  the  interaction  of  the  sharp  change  in  topography  and  the  strong 
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tidal  currents  near  the  Luzon  Strait.  The  largest  ever  observed  internal  waves 
occur  in  this  region  (Hsu  &  Liu  2000).  These  waves  dominate  the  ocean 
variability  in  the  northeastern  (NE)  SCS  basin.  These  non-linear  waves  can 
significantly  effect  the  acoustic  propagation  in  the  basin. 

C.  WISE  EXPERIMENT 

The  Windy  Islands  Soliton  Experiment  (WISE)  was  an  integrated 
acoustics  and  physical  oceanography  experiment,  an  international  collaboration 
between  the  Naval  Postgraduate  School  (NPS),  National  Taiwan  University, 
National  Sun  Yet-Sen  University,  Ocean  Acoustical  Services  and  Instrumentation 
Systems  (OASIS),  and  Woods  Hole  Oceanographic  Institute  (WHOI).  The 
experiment  took  place  in  the  NE  SCS  and  Luzon  Strait  and  lasted  one  year,  April 
2005-2006.  One  of  the  main  goals  of  WISE  was  to  study  the  evolution  of  trans¬ 
basin  nonlinear  internal  waves  (Ramp  et  al.  2010).  Moorings  were  placed  on  the 
shelf,  in  the  Luzon  strait,  and  the  deep  basin. 

In  WISE  two  acoustic  propagation  studies  were  conducted,  one  on  the 
continental  shelf  and  one  in  the  deep  basin.  The  focus  of  these  was  to  study  the 
physics,  phenomenology,  and  statistics  of  sound  propagation  in  the  NE  SCS. 
The  objective  of  the  basin  propagation  experiment  was  to  study  and  characterize 
the  supertidal  to  seasonal  scale  impacts  of  the  trans-basin,  nonlinear  internal 
waves  on  long-range  transmission  loss,  and  to  help  monitor  the  evolution  of  the 
trans-basin  internal  waves  in  the  basin’s  interior  (Miller  et  al.  2009).  This  was 
completed  by  the  use  of  two  deep-water  moorings  that  transmitted  and  received 
a  400  Hz  signal  repeatedly  over  a  one-year  period.  Modeling  of  the  observed 
arrival  structure  is  the  focus  of  this  thesis.  The  location  of  these  moorings  are 
shown  in  Figure  1 . 
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Figure  1 .  Locations  of  the  source  (B2)  and  receiver  (B1 )  in  the  SCS 

superimposed  on  the  bathymetry  (From  Bernotavicius  et  al.  2010). 

D.  OBJECTIVE 

The  objective  of  this  research  is  to  model  and  study  the  effects  of  the 
transbasin  internal  waves  on  low-frequency  signal  transmission  through  the  SCS 
basin.  Specifically,  the  fluctuations  in  the  multipath  arrival  structure  of  a  400-Hz 
acoustic  pulse  transmitted  through  a  distance  of  167-km  in  the  NE  SCS  basin  in 
the  presence  an  internal  ocean  soliton.  Previous  research  by  Bernotavicius 
(2010)  modeled  the  background  arrival  structure  in  the  unperturbed  SCS  basin 
using  the  Hamiltonian  Acoustic  Ray-tracing  Program  for  the  Ocean  (HARPO). 
That  research  used  a  silty  clay  type  bottom.  A  bottom  attenuation  rate  of  .486 
dB/km/Hz,  a  bottom  density  of  1.24  kg/km3  and  a  bottom  sound  speed  of  1,521 
m/s  were  used.  That  accurately  modeled  the  mean  arrival  structure  for  the  SCS 
basin.  An  initial  launch  angle  angular  resolution  of  .005  degrees  was  used.  That 
resolution  or  finer  was  able  to  capture  all  the  eigenrays  (Bernotavicius  et  al. 
2010).  This  research  will  study  the  fluctuations  in  the  arrival  structure  caused  by 
the  nonlinear  ocean  solitons. 
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E.  APPROACH 

The  approach  of  this  thesis  was  to  develop  a  numerical  model  of  a  400  Hz 
signal  transmission  through  the  SCS  in  the  presence  of  a  trans-basin  internal 
soliton  and  then  examine  and  contrast  the  modeled  arrival  structures  at  the 
receiver.  The  model  is  based  on  ray  theory,  and  entails  the  integration  of  a 
raytracing  model  and  an  eigenray  ray  search  and  arrival  structure  simulation 
code.  The  model  permits  the  study  of  the  physics  and  phenomenology  of  sound 
propagation  though  the  SCS  trans-basin  nonlinear  internal  waves.  A  modeled 
range  dependent  sound  speed  profile,  measured  bathymetry  and  inferred 
bottom-loss  characteristics  from  previous  research  were  used. 

F.  THESIS  OUTLINE 

The  rest  of  this  thesis  contains  three  chapters.  Chapter  II  describes  the 
ray  trace  model,  the  eigenray  search  and  arrival  structure  simulation  model,  and 
the  methodology  used  to  integrate  the  ray  trace  model  with  the  eigenray  search 
arrival  structure  simulation.  Chapter  III  describes  the  results  of  the  five  cases 
studied.  A  summary  of  results  and  conclusions  are  given  in  Chapter  IV. 
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II.  METHODOLOGY 


A.  ACOUSTIC  MODEL 

1 .  Ray  T race  Theory 

Modeling  of  sound  propagation  through  the  ocean  requires  solving  the 


governing  spherical  wave  equation: 


V2p- 


1  d2p 


dt2 


p) 


■  0 


(2.1) 


subject  to  the  appropriate  surface  and  bottom  boundary  conditions.  Here  p  is  the 
sound  pressure,  c  is  sound  speed,  t  is  time,  and  r  is  a  position  vector  (Medwin, 
2005).  The  assumption  is  made  in  ray  acoustics  that  variation  in  sound  speed  is 
gradual  over  a  small  number  of  acoustic  wavelengths.  It  is  also  assumes  that 
solutions  of  the  equation  take  the  following  form: 


p(r,t)  =  A(r)exp^2nft-k0W(r)  j 


(2.2) 


where  L 


2  nf 


is  the  reference  wave  number,  c0  is  the  reference  sound  speed, 


/  is  the  acoustic  frequency,  and  W (V)  is  the  eikonal  when  multiplied  by  k0  gives 

the  phase  component  of  the  solution.  Equation  2.2  is  substituted  into  Equation 
2.1  to  yield  the  eikonal  and  transport  equations  in  Equation  2.3.  The  eikonal 
yields  the  ray  equations  and  the  ray  paths  while  the  transport  equation  governs 
the  ray  amplitude.  The  phase  is  determined  by  integrating  ekional  along  the  ray 
path.  This  results  in  the  phase  equaling  2/t/t  where  t  is  the  ray  travel  time. 


The  ratio  of  a  reference  sound  speed  to  sound  speed 


'(r) 


is  the  quantity  n, 


called  the  index  of  refraction  and  s  is  the  distance  along  the  ray  path. 
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Eikonal  Equation:  |VlY(r)|  =  n2(f ) 

Transport  Equation:  2VA(r)- WW (r)  +  A(r)V2W (r)  =  0  (2.3) 

\ 

=  Vn(r) 

Ray  theory  assumes  that  diffractive  leakage  is  small  and  the  change  in 
amplitude  occurs  slowly.  The  transport  equation  implies  that  power  must  be 
conserved  along  a  ray  tube.  A  full  derivation  of  these  equations  can  be  found  in 
a  reference  such  as  Ziomek  (1995). 

A  three-dimensional  ray  trace  model  developed  by  Y.  T.  Lin  (2012, 
personal  communication)  at  WHOI  was  used  for  this  thesis.  This  model  is 
MATLAB  based  and  is  intended  to  replace  the  HARPO  model  used  in  previous 
research.  HARPO  is  a  FORTRAN  based  model  that  can  only  run  on  obsolete 
computer  systems.  The  model  was  reduced  to  a  two  dimensional  model.  An 
initial  launch  angle  angular  resolution  of  .005  degrees  was  used  in  the  raytracing. 
An  initial  launch  angle  angular  aperture  of  -20  to  20  degrees  was  used.  This 
captures  all  eigenrays  that  significantly  contribute  and  minimizes  computational 
requirements  (Bernotavicius  et  al.  2010). 

2.  Eigenray  Determination 

To  model  the  multipath  acoustic  arrival  structure,  it  is  necessary  to  identify 
the  acoustic  paths  connecting  the  source  to  the  receiver,  i.e.,  eigenrays.  The 
eigenray  search  and  arrival  structure  simulation  MATLAB  code  of  Chiu  (1994) 
was  integrated  with  output  of  the  raytrace  model  to  calculate  the  eigenray 
parameters  and  arrival  structure.  This  code  determines  the  depth  difference  at 
the  receiver  range  between  the  receiver  and  the  ray  as  a  function  of  initial  launch 
angle.  Locating  the  zeros  on  this  curve  allows  for  the  eigenray  parameters  to  be 
interpolated  accurately.  No  rays  are  retraced  while  determining  the  eigenrays. 
For  the  plotting  of  eigenrays,  those  rays  located  nearest  to  the  eigenrays  are 
used.  This  code  accurately  interpolates  using  higher  order  polynomial  fitting  for 
the  eigenray  parameters  including  amplitude,  phase,  and  phase  front  arrival  time. 


Ray  Equation :  —  n[r)  — 
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Once  the  eigenray  parameters  are  determined  the  arrival  structure  can  be 
generated  using  Equation  2.4.  A  model  pulse  is  created  to  match  the  energy  and 
shape  of  the  source  signal  complex  envelop  used  in  WISE.  Given  this  pulse 
complex  envelop  and  the  eigenray  parameters  the  arrival  structure  can  be 
calculated  as: 

?{<)=  I  (2.4) 

eigenrays 

where  f  represents  the  complex  envelop  at  the  receiver  and  s  represents  the 
complex  envelop  at  the  source.  A  more  detailed  discussion  of  this  approximation 
can  be  found  in  Hager  (2008). 

B.  MODEL  INPUT 

1.  Model  Domain 

In  the  model,  the  sound  source  and  receiver  were  located  at  the  B2  and 
B1  locations  shown  in  Figure  1,  placed  167  km  apart  and  located  at  a  depth  of 
800  m.  This  was  done  to  model  the  sound  source  at  the  B2  mooring  and  the 
receiver  at  the  B1  mooring.  A  signal  center  frequency  of  400  Hz  was  used.  A 
bottom  attenuation  rate  of  .486  dB/km/Hz,  a  bottom  density  of  1 .24  kg/km3,  and  a 
sediment  sound  speed  of  1 .521  km/s  were  used.  These  bottom  parameters  best 
represent  the  bottom  loss  characteristics  of  the  NE  SCS  (Bernotavicius  et  al. 
2010).  An  rms  surface  roughness  of  zero  was  chosen  to  represent  a  smooth  sea 
surface.  This  was  done  to  remove  any  losses  due  to  surface  roughness  suffered 
by  the  rays.  Bathymetry  was  obtained  from  high  resolution  transects  of  the 
Luzon  Strait  (Liu  et  al.  1998).  This  data  was  smoothed  using  a  cubic  spline  to 
give  a  continuous  bathymetry  for  raytracing. 

2.  Sound  Speed 

For  this  research  a  two  dimensional  range  dependent  sound  speed  profile 
was  used.  The  modeled  perturbed  sound  speed  profile  consists  of  a  background 
unperturbed  profile  c(z)  that  is  range  independent  and  a  perturbation  8c(r,z), 
which  depends  on  range  r  and  depth  z. 
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c{r,z)  =  c(z)  +  8c{r,z)  (2.5) 

The  sound  speed  profile  was  determined  using  the  Medwin  sound  speed 
equation  (Medwin,  1975)  applying  temperature  salinity  and  depth  data.  The  data 
includes  a  three-month  time  series  of  temperature  with  a  sampling  rate  of  one 
minute,  at  the  B1  mooring.  This  mooring  provided  temperatures  at  the  following 
nine  depths  in  meters:  125,  173,  375,  500,  755,  1000,  1385,  1884,  and  2438. 
The  data  also  includes  the  mean  salinity  profile  from  the  CTD  casts  conducted 
during  the  WISE  experiment.  The  minimum  value  of  sound  speed  for  each  depth 
was  used  to  represent  the  unperturbed  sound  speed  profile.  This  is  due  to  the 
fact  that  the  transbasin  solitons  are  depression  internal  waves  and  hence  they 
would  only  lead  to  sound  speed  increases.  For  depths  less  than  125  meters,  the 
sound  speed  gradient  at  125  meters  was  extrapolated  to  the  surface.  For  depths 
below  2438  meters,  the  sound  speed  was  extrapolated  using  a  constant  sound 
speed  gradient  of  .016  s"1 . 

3.  Sound  Speed  Perturbation 

The  sound  speed  perturbations  were  assumed  to  have  a  separable  form, 
consisting  of  a  vertical  component  and  a  horizontal  component.  The  equation  for 
the  perturbation  is  therefore. 

8c{r,z)  =  a(z)’  fi(r)  (2.6) 

To  determine  the  vertical  component  of  the  sound  speed  perturbations,  an 
Empirical  Orthogonal  Function  (EOF)  analysis  was  performed  on  the  sound 
speed  time  series  calculated  from  the  temperature,  salinity  and  depth 
information.  EOF  analysis  is  a  common  tool  used  to  determine  the  dominant 
structure  of  a  data  set.  EOF  analysis  decomposes  the  data  into  orthogonal 
functions.  The  resulting  orthogonal  functions  will  be  rank  ordered  by  variance. 
The  first  orthogonal  function  contains  the  most  variance.  The  EOF  analysis 
conducted  on  the  sound  speed  time  series  showed  that  the  first  orthogonal 
function  contained  95  percent  of  the  variance.  The  remaining  five  percent  of  the 
variance  was  spread  among  the  remaining  eight  orthogonal  functions.  The  first 
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orthogonal  function  was  chosen  to  represent  the  vertical  structure  of  a  soliton.  A 
more  detailed  description  of  this  method  of  data  analysis  can  be  found  in  Wilks 
(2006).  The  vertical  component  is  shown  in  Figure  2.  This  was  obtained  by  a 
cubic  spline  fit  through  the  points  of  the  first  orthogonal  function  and  zero  at  the 
surface  and  a  depth  of  4,000  m.  The  vertical  structure  was  normalized  to  have  a 
maximum  value  of  unity.  The  depth  of  maximum  influence  for  the  modeled 
soliton  is  91  m. 
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Figure  2.  Sound  speed  perturbation  vertical  structure  derived  from  EOF  analysis 

of  sound  speed  profile  time  series. 


The  modeled  horizontal  component  is  shown  in  the  Equation  2.7,  which  is 
adopted  from  (Chiu  et  al.  2004): 


[3(r)  =  8  sech2 


f 

In 

v 


r~ro' 

A 


8  =  15  — 

s 

A  =  11,520  m 
r0  =  Center  location 


(2.7) 
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where  8  is  the  soliton  amplitude,  A  is  the  soliton  width,  and  r0  is  the  soliton 
center  location.  The  values  used  for  the  model  soliton  width  and  amplitude  were 
taken  from  the  WISE  dataset.  The  soliton  width  was  shifted  from  time  to  space 
using  a  soliton  phase  speed  of  1.25  m/s  (Ramp  et  al.  2004).  How  closely  this 
structure  matches  observed  data  is  shown  in  Figure  3.  For  this  research,  a 
soliton  was  placed  in  five  different  locations  40,  60,  80,  120,  and  160  km  from  the 
source.  This  is  to  represent  a  soliton  transiting  the  basin.  The  vertical  and 
horizontal  components  are  combined  to  represent  the  sound  speed  perturbation 
that  is  shown  in  Figure  4. 


Figure  3.  The  horizontal  structure  of  the  sound  speed  perturbation.  Both  the 

modeled  (red)  and  observed  (blue)  are  shown. 
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Case  1  Case  2  Case  3 
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Figure  4.  Contour  plot  of  the  modeled  sound  speed  with  an  internal  soliton 
located  at  40  km  overlaying  the  bathymetry.  The  bottom  attenuation 
rate,  bottom  density  and  bottom  sound  speed  are  also  displayed. 


For  this  research  a  soliton  was  placed  at  five  different  locations  in  the 
model  domain.  A  model  run  was  also  conducted  with  no  soliton;  this  will  be 
referred  to  as  the  unperturbed  case.  The  collection  of  the  five  cases  is  meant  to 
represent  a  soliton  transiting  across  the  NE  SCS  basin.  The  location  of  the 
soliton  center  in  each  case  is  shown  in  Table  1 . 


Case  # 

Location  of  Soliton  Center 

1 

40  km 

2 

60  km 

3 

80  km 

4 

120  km 

5 

160  km 

Table  1 .  Center  location  of  the  ocean  internal  soliton  for  each  case. 
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III.  RESULTS 


A.  OVERVIEW 

Consistently,  each  perturbation  case  and  the  unperturbed  case  has  six 
separate  arrivals.  Their  arrival  times  are  approximately  constant  at:  112.8s, 
112.9s,  113.1s,  113.3s,  113.5s,  and  113.7s,  respectively.  A  plot  of  the  modeled 
sound  pressure  (in  dB  re  1  /uPa)  versus  travel  time  for  each  of  the  five  cases  and 
the  unperturbed  case  is  shown  in  Figure  5.  Although  the  arrival  times  are  rather 
robust,  large  fluctuations  in  amplitude  for  each  arrival  are  present.  The  factors 
that  affect  the  amplitude  for  each  of  the  arrivals  are:  bottom  interaction,  change  in 
ray  tube  spreading,  phase  interference  and  an  increase/decrease  in  the  number 
of  eigenrays  making  up  an  arrival.  These  factors  will  be  illustrated  in  the 
subsequent  sections.  In  two  of  the  cases  the  second  arrival  diminishes  greatly. 
Each  arrival  for  each  perturbation  case  is  slightly  ahead  of  the  corresponding 
arrival  in  the  unperturbed  case.  This  is  due  to  the  transbasin  soliton  being  a 
depression  wave  that  leads  to  an  increase  in  sound  speed. 


—Unperturbed 


Figure  5.  Acoustic  arrival  structure  (receive  level  versus  travel  time)  for  each 

case  and  the  unperturbed  case. 
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For  some  arrivals  in  the  perturbed  cases  the  associated  number  of 
eigenrays  changes.  This  is  associated  with  micropathing  caused  by  the  relatively 
strong  gradients.  In  addition  the  sharp  sound  speed  gradient  in  the  perturbed 
cases,  the  locations  of  bottom  interactions  can  deviate  slightly  and  thus  the  same 
group  of  rays  can  experience  different  bathymetry  gradients.  An  example  of  the 
existence  of  micropathing  can  be  seen  in  Figure  6.  Where  the  curve  of  arrival 
depth  difference  (i.e.,  arrival  depth  at  the  receiver  range  minus  the  receiver 
depth)  versus  the  initial  launch  angle  of  rays  is  smooth,  near  a  low  initial  launch 
angle,  no  micropathing  exists.  As  the  initial  launch  angle  increases  and  the 
depth  difference  versus  initial  launch  angle  curve  sharply  changes  micropathing 
begins  to  develop.  Note  that  the  zero  crossings  of  this  curve  define  the  initial 
launch  angle  of  all  the  eigenrays.  The  remainder  of  this  chapter  will  explore  the 
physics  and  phenomenology  of  the  amplitude  changes  for  each  case. 


2 

1.5 

1 


"^5  -10  -5  0  5  10  15 

Initial  Launch  Angle  (deg) 


Figure  6.  Arrival  depth  versus  initial  launch  angle.  The  depth  difference  is  the 
deviation  from  the  receiver  depth.  The  zeros  crossings  are  the  initial 
launch  angles  of  the  eigenrays  (shown  as  green  stars). 


14 


B.  CASE  1 


For  Case  1  the  soliton  is  placed  at  40  km  from  the  source.  A  plot  of  the 
acoustic  arrival  structure  for  this  case  is  shown  in  Figure  7.  The  first  and  fourth 
arrivals  have  a  sound  pressure  lower  than  the  unperturbed  case.  The  second 
arrival  for  case  1  has  a  higher  level  than  the  unperturbed  case.  The  third,  fifth, 
and  sixth  arrivals  are  roughly  the  same. 


120 


-Unperturbed 
-40  KM  Perturbation 


114 


Figure  7.  Acoustic  arrival  structure  (receive  level  versus  travel  time)  for  Case  1 

(blue)  and  for  the  unperturbed  case  (red). 

The  ray  paths  of  the  eigenrays  for  the  first,  second  and  fourth  arrivals  are 
shown  in  Figure  8.  Case  1  is  in  blue,  and  the  unperturbed  case  is  in  red.  In  the 
first  arrival,  a  shift  in  the  initial  launch  angles  is  seen.  For  the  second  arrival,  the 
soliton  causes  more  eigenrays  to  form,  leading  to  a  higher  receive  level.  Case  1 
has  four  eigenrays  while  the  unperturbed  case  only  has  two.  The  fourth  arrival 
shows  that  the  eigenrays  for  the  perturbed  case  have  an  additional  bottom 
interaction  located  around  20  kilometers  from  the  source  that  the  unperturbed 
case  does  not  show.  This  is  highlighted  with  a  black  circle  in  Figure  8. 
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Figure  9  is  a  plot  of  ray  tube  spreading  loss  versus  initial  launch  angle  for 
Case  1  and  the  unperturbed  case.  The  discrete  symbols  represent  the 
eigenrays.  A  shift  in  the  structure  of  this  plot  can  clearly  be  seen.  The  first 
arrival  has  increased  raytube  spreading  loss  for  several  of  the  eigenrays  in  Case 
1 .  This  is  highlighted  with  a  black  circle  in  Figure  9. 


Figure  8.  The  eigenrays  for  the  1st,  2nd,  and  4th  arrivals  for  Case  1  (blue)  and  for 
the  unperturbed  case  (red).  Additional  bottom  interactions  for  4th 
arrival  of  Case  1  are  highlighted  in  the  black  circle. 
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Figure  9.  Raytube  spreading  loss  versus  initial  launch  angle  (solid  lines)  for 
Case  1  (blue)  and  for  the  unperturbed  case  (red).  Discrete  symbols 
denote  the  eigenrays.  Increased  spreading  loss  for  the  1st  arrival  of 
Case  1  is  highlighted  in  the  black  circles. 


A  plot  of  the  number  of  bottom  bounces  versus  initial  launch  angle  is 
shown  in  Figure  10.  Case  1  is  in  blue  and  the  unperturbed  case  is  in  red.  The 
fourth  arrival  for  Case  1  contains  several  eigenrays  with  5  bottom  bounces  while 
the  same  arrival  for  unperturbed  case  has  only  4  bottom  bounces.  This  is  shown 
with  a  black  circle  in  Figure  10.  This  increase  in  bottom  interactions  leads  to  an 
increase  bottom  loss  for  the  fourth  arrival  seen  in  Figure  1 1 . 
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Figure  10.  Number  of  bottom  bounces  versus  initial  launch  angle  for  Case  1 

(blue)  and  for  the  unperturbed  case  (red).  Additional  bottom  bounces 
for  the  4th  arrival  of  Casel  are  highlighted  with  black  circles. 
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Figure  1 1 .  Bottom  loss  versus  initial  launch  angle  (solid  lines)  for  Case  1  (blue) 
and  for  the  unperturbed  case  (red).  Discrete  symbols  denote 

eigenrays. 


18 
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CASE  2 


For  Case  2,  the  soliton  is  placed  at  60  km  from  the  source.  A  plot  of  the 
acoustic  arrival  structure  for  this  case  is  shown  in  Figure  12.  The  second  and 
sixth  arrivals  are  roughly  10  dB  lower  than  the  unperturbed  case.  The  first,  third, 
fourth,  and  fifth  arrivals  are  all  similar  to  the  unperturbed  case. 
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Figure  12.  Acoustic  arrival  structure  (receive  level  versus  travel  time)  for  Case  2 

(blue)  and  for  the  unperturbed  case  (red). 


The  eigenrays  for  the  second  and  sixth  arrivals  for  Case  2  are  shown  in 
Figure  13.  Case  2  is  in  blue  and  the  unperturbed  case  in  red.  The  second  arrival 
for  Case  2  and  the  unperturbed  case  have  only  two  eigenrays.  These  two  rays 
for  Case  2  are  more  out  of  phase  then  the  rays  of  unperturbed  case  and 
therefore  have  destructive  interference  for  this  arrival.  The  sixth  arrival  for  Case 
2  consists  of  13  eigenrays  while  the  same  arrival  for  the  unperturbed  case  has  15 
eigenrays.  This  difference  in  the  number  of  eigenrays  combined  with  the  phase 
interference  explains  the  lower  sound  pressure. 
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Figure  13.  The  eigenrays  for  the  2nd  and  6th  arrivals  for  Case  2  (blue)  and  for  the 

unperturbed  case  (red). 


A  plot  of  ray  tube  spreading  loss  versus  initial  launch  angle  for  Case  3  and 
for  the  unperturbed  case  is  shown  in  Figure  14.  A  plot  of  the  corresponding 
bottom  loss  versus  initial  launch  angle  is  shown  in  Figure  15.  The  second  arrival 
for  Case  3  and  the  unperturbed  case  have  similar  values  for  both  ray  tube 
spreading  loss  and  bottom  loss.  This  is  shown  with  blacks  circles  in  Figure  14 
and  15.  This  indicates  that  destructive  interference  caused  by  the  two  rays  being 
more  out  of  phase  is  the  cause  of  the  decreased  sound  pressure  for  the  second 
arrival  of  Case  3. 
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Figure  14.  Raytube  spreading  loss  versus  initial  launch  angle  (solid  lines)  for 
Case  2  (blue)  and  for  the  unperturbed  case  (red).  Discrete  symbols 
denote  eigenrays.  Black  circles  highlight  similar  values  for  2nd  arrival 
for  Case  2  and  the  unperturbed  case. 
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Figure  15.  Bottom  loss  versus  initial  launch  angle  (solid  lines)  for  Case  2  (blue) 
and  for  the  unperturbed  case  (red).  Discrete  symbols  denote 
eigenrays.  Black  circles  highlight  similar  values  for  2nd  arrival  for  Case 

2  and  the  unperturbed  case. 
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D.  CASE  3 

For  Case  3,  the  soliton  is  placed  at  80  km  from  the  source.  A  plot  of  the 
acoustic  arrival  structure  for  this  case  is  shown  in  Figure  16.  The  first,  second, 
and  sixth  arrivals  have  a  sound  pressure  than  the  unperturbed  case.  The  fourth 
arrival  for  case  3  has  a  higher  sound  pressure  than  the  unperturbed  case. 
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Figure  16.  Acoustic  arrival  structure  (receive  level  versus  travel  time)  for  case  3 

(blue)  and  for  the  unperturbed  case  (red). 

The  eigenrays  for  the  first,  second,  fourth,  and  sixth  arrivals  for  Case  3  are 
shown  in  Figure  17.  Case  3  is  in  blue  and  the  unperturbed  case  is  in  red.  The 
first  arrival  is  lower  due  to  that  the  eigenrays  have  additional  bottom  interactions. 
For  the  second  arrival  there  are  two  eigenrays  for  Case  3  and  the  unperturbed 
case.  The  two  rays  are  more  out  of  phase  for  Case  3  than  the  unperturbed  case. 
The  sixth  arrival  for  case  3  contains  only  8  eigenrays,  while  the  unperturbed  case 
contains  15  rays. 
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Figure  17.  The  eigenrays  for  the  1st,  2nd,  4th,  and  6th  arrivals  for  Case  3  (blue)  and 

for  the  unperturbed  case  (red). 


A  plot  of  the  number  of  bottom  bounces  versus  initial  launch  angle  is 
shown  in  Figure  18.  Several  of  the  eigenrays  for  the  first  arrival  of  Case  3  have 
additional  bottom  interactions.  This  increased  bottom  interaction  lead  to 
increased  bottom  loss  for  the  first  arrival  of  Case  3.  A  plot  of  bottom  loss  versus 
initial  launch  angle  is  shown  in  Figure  19.  The  increased  bottom  loss  is 
highlighted  with  a  black  circle  in  Figure  19. 
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Figure  18.  Number  of  bottom  bounces  versus  initial  launch  angle  for  Case  3 

(blue)  and  for  the  unperturbed  case  (red). 
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Figure  19.  Bottom  loss  versus  initial  launch  angle  (solid  lines)  for  Case  3  (blue) 
and  for  the  unperturbed  case  (red).  Discrete  symbols  denote 
eigenrays.  Black  circles  highlight  the  higher  bottom  loss  for  the  1st 

arrival  of  Case  3. 
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The  fourth  arrival  for  Case  3  has  less  bottom  loss  than  the  unperturbed 
case.  This  is  highlighted  with  a  black  circle  in  Figure  20.  This  change  is  due  to 
the  change  of  the  ray  angle  at  bottom  for  these  rays.  The  eigenrays  for  the 
fourth  arrival  are  represented  as  triangles. 
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Figure  20.  Ray  tube-spreading  loss  versus  initial  launch  angle  (solid  lines)  for 
Case  3  (blue)  and  for  the  unperturbed  case  (red).  Discrete  symbols 
denote  eigenrays.  Black  circle  highlight  the  difference  in  spreading 

loss  for  the  4th  arrival. 


E.  CASE  4 

For  case  4,  the  introduction  of  a  soliton  did  not  significantly  change  the 
arrival  structure.  In  case  4,  the  soliton  is  placed  at  120  km  from  the  source.  A 
plot  of  the  acoustic  arrival  structure  for  this  case  is  shown  in  Figure  21.  The 
eigenrays  for  case  4  follow  a  very  similar  path  to  the  eigenrays  of  the 
unperturbed  case.  A  plot  of  the  eigenrays  for  all  of  the  arrivals  in  case  4  is  shown 
in  Figure  22. 
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Figure  21 .  Acoustic  arrival  structure  (receive  level  versus  travel  time)  for  Case  4 

(blue)  and  for  the  unperturbed  case  (red). 
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Figure  22.  The  eigenrays  for  all  arrivals  for  Case  4  (blue)  and  for  unperturbed 

case  (red). 
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F.  CASE  5 

For  Case  5,  the  soliton  is  placed  at  160  km  from  the  source.  A  plot  of  the 
acoustic  arrival  structure  for  this  case  is  shown  in  Figure  23.  The  arrival  structure 
for  Case  5  is  very  similar  to  the  unperturbed  case.  The  soliton  has  very  little 
effect  on  the  eigenrays.  A  plot  of  the  eigenrays  for  all  of  the  arrivals  for  Case  5  is 
shown  in  Figure  24. 


Figure  23.  Acoustic  arrival  structure  (receive  level  versus  travel  time)  for  Case  5 

(blue)  and  for  the  unperturbed  case  (red). 
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Figure  24.  The  eigenrays  for  all  arrivals  for  Case  5  (blue)  and  for  the  unperturbed 

case  (red). 
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IV.  CONCLUSIONS 


The  objective  of  this  research  was  to  model  and  study  the  effects  of  the 
transbasin  internal  waves  on  low-frequency  signal  transmission  through  the  SCS 
basin.  Specifically,  the  fluctuations  in  the  multipath  arrival  structure  of  a  400-Hz 
acoustic  pulse  transmitted  through  a  distance  of  167-km  in  the  NE  SCS  basin  in 
the  presence  an  internal  ocean  soliton  were  studied.  The  approach  of  this  thesis 
was  to  develop  a  numerical  model  of  a  400  Hz  signal  transmission  through  the 
SCS  in  the  presence  of  a  trans-basin  internal  soliton  and  then  examine  and 
contrast  the  modeled  arrival  structures  at  the  receiver.  The  model  is  based  on  ray 
theory,  and  entails  the  integration  of  a  raytracing  model  and  an  eigenray  ray 
search  and  arrival  structure  simulation  code.  The  model  permits  the  study  of  the 
physics  and  phenomenology  of  sound  propagation  though  the  SCS  trans-basin 
nonlinear  internal  waves.  A  modeled  range  dependent  sound  speed  profile, 
measured  bathymetry  and  inferred  bottom-loss  characteristics  from  previous 
research  were  used. 

Each  perturbation  case  and  the  unperturbed  case  have  six  separate 
arrivals  with  robust  travel  time.  However,  large  amplitude  fluctuations  were 
observed  for  each  of  these  arrivals.  The  factors  that  affect  the  amplitude  for 
each  of  the  arrivals  are  changes  in:  the  number  of  bottom  interactions,  ray  tube 
spreading,  phase  interference  and  in  the  number  of  eigenrays  making  up  an 
arrival.  Micro-multipath  arrivals  affect  a  number  of  eigenrays  for  all  arrivals.  The 
results  also  show  that  the  closer  the  soliton  to  the  receiver,  the  less  impact  the 
soliton  had  on  the  arrival  structure.  This  is  because  when  the  soliton  was  placed 
farther  from  the  receiver,  more  bottom  interactions  downstream  existed. 
Therefore  changing  the  locations  of  these  bottom  interactions  resulted  in  more 
changes  in  the  arrival  amplitude.  As  the  range  between  the  receiver  and  the 
soliton  was  decreased  the  perturbation  had  less  effect  on  the  location  of  bottom 
interactions  and  cause  less  change  in  the  arrival  structure. 
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Future  research  will  entail  a  detailed  model-data  comparison  of  this  model 
and  the  observations  from  the  WISE  experiment.  This  comparison  will  be  both 
qualitative  and  quantitative  in  nature  to  further  validate  the  model. 
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